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Abstract. 

The survival of qq bound states at temperatures higher than the chiral restoration temperature, T c , 
recently observed in lattice QCD, is discussed in the framework of the Nambu Jona-Lasinio model. 
The perturbative determination of the spectral function provides an indication of a pseudogap phase 
above T c . 
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IT) ■ Recent lattice studies about the Quantum Chromodynamics transition at finite tempera- 

ture from the hadronic to the deconfined phase, show that, heavy, and possibly also light 
mesonic bound states survive up to twice the deconfinement temperature (for a recent re- 
view see Jl|] and references therein), in contrast with the former suggestion that a strong 
suppression of heavy bound states, such as the J /\ff, would occur just above the critical 
temperature. This unexpected feature indicates some not yet understood mechanism at 
the deconfinement transition. In addition, the first experimental results collected at RHIC 
support the presence of a strongly interacting phase above the critical temperature. 

This picture has some analogies with the physics of high temperature superconductors 
where the coherence length is much smaller than in ordinary superconductors and the 
temperature corresponding to the onset of superconductivity turns out to be much lower 
than the temperature related to the Cooper pairs formation. Between these two tempera- 
tures a pseudogap is observed, which consists in a depletion of the single particle density 
of states around the Fermi level. These features could be regarded as a phenomenologi- 
cal manifestation of a crossover from the ordinary Bardeen-Cooper-Schrieffer supercon- 
ducting behavior to a Bose-Einstein Condensate behavior. 

An explanation of the high temperature superconductivity has been provided by 
Emery and Kivelson (see [2]) who suggested that phase fluctuations of the condensate 
should be responsible for the spoiling of the long range coherence, the appearance 
of the pseudogap and the consequent lowering of the temperature associated to the 
onset of superconductivity . The relation between the pseudogap transition and field 
fluctuation has been explored within the framework of the Nambu Jona-Lasinio (NJL) 
model Oiyiy], and also indicated as a possible explanation of the observed lattice bound 
states above the critical temperature !□]• 



Following |a], after a brief analysis of some results in mean field theory for the NJL 
model, we briefly discuss the effects of the fluctuations and the consequent appearance of 
two characteristic temperatures: T c , corresponding to the restoration of chiral symmetry 
for the NJL model, and T* associated to the decoupling of the mesons qq from the 
spectrum. Finally we present a calculation of the pseudogap between T c and T*. 



NJL MODEL IN THE MEAN FIELD APPROXIMATION 



We consider the NJL model with an isospin doublet of massless quarks with three 
colors (Nf = 2, N c = 3 ) at finite chemical potential ji 
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As it is well known (for reviews on the NJL model see e.g. UilO]), the mean field 
approximation provides the self-consistent equation (at jl = 0) for the mass gap m q : 
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where A is a 3D cutoff. Moreover the pion decay constant f n , again at /i = 0, is given 
by 

fl = -Aim 2 a N c [ -^-r—i \ ^ (3) 

and one can use Eqs. © and © to get A by fixing m q . In particular we use m q — 300 
MeV and f K = 93.3 MeV as input, which gives A = 675 MeV. 

The approach is then easily generalized at finite temperature and density to obtain a 
new gap equation for m q (T,ji) which, after getting rid of the coupling Go by means of 
Eq. ©, reads 
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where e = J p 2 + m q (J ', jii) and x = }i/T and y = e/T. 

The gap equation at finite temperature and density naturally provides a critical line 
signaling a phase transition. In fact, the vanishing of the solution of Eq. ©: m q (T, fl) = 
defines the critical temperature T*(fi) and its plot in the T — jl plane is displayed in fig. 
1, with the input for the 3D cutoff A obtained from the values of m q and f n at T = fl = 
given above. Our next step is to check to what extent the fluctuations can modify the 
results of the mean field analysis. 
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FIGURE 1. Plot of T* (ju ) in MeV, as obtained from Eq. (gj for A = 675 Me V. 



FLUCTUATION EFFECTS 

In order to study the modifications to the mean field analysis it is necessary to introduce 
a method that can properly account for the fluctuation effects. Such approach has been 
developed 0,01 and its essential features will be briefly illustrated in the following. 

The four fermion interaction in the model in Eq. ([T]) can be cast in the form of Yukawa 
coupling, by introducing the auxiliary scalar and pseudoscalar fields a and % : 

J2f = Wdvf + ny -go(o + iy 5 T-7i))Y-^ [ a 2 + jp] (5) 

with the quark-meson coupling constant given by go = m q /f K , which is the analogous 
of the Goldberger-Treiman relation. By functional integration of the fermionic degrees 
of freedom and derivative expansion, Eq. © becomes equivalent to the o model with 
lagrangian 



^ = § ((^) 2 + [dnf -^(o 2 + ie-flf) 



(6) 



where the stiffness parameter, /3, which is relevant for our purposes, is given by 
j m \ j f / dA PE 1 1 f d 4 p E p 2 E 



P - WcNf [ ^ {2n)4 {p 2 +m 2 )2 2 J A (2^)4 (pi +m2)3 ) (?) 

and pe is the Euclidean four-momentum. 

The key point is that an equivalent description of the problem should be obtained by 
resorting to an effective theory, defined by a non-linear o model for the fields o and % , 
satisfying the constraint o 2 + 7Z 2 = f% [4, 5]. This constraint can be introduced into the 
functional generator Z of the non-linear a model, by means of a functional Lagrange 
multiplier A (x) 

Z = J [do] [dn] [dX] exp j j| j d A x ((do) 2 + (dn) 2 - X (o 2 + % 2 - f 2 ) ) j (8) 



The parameter /3 in Eq. © indicates the stiffness of the non-linear o model. 

The functional integration of the % fields in Eq. ® can be performed and, afterward, 
one can use the saddle point approximation and search for x-independent solutions for 
o and A. As it is evident from Eq. ©, the constant A plays the role of a square mass 
both for o and % fields. Since we are looking for phase transitions at finite T and /!, 
we consider Matsubara frequencies cD„ = lltnT and non-vanishing baryonic chemical 
potential and the two saddle point conditions read 



= Aa; = P(f 2 -a 2 )-(N 2 -l) £ / 



<Pk 



J±Ja* {Ik? (a)n-i^) 2 +k 2 + X 



(9) 



where is a suitable cutoff for the non-linear o model. 

The first condition = Act indicates that at least one of the two variables must vanish 
and the other variable is implicitly determined from the second condition in Eq. ©, in 
terms of T, /i and of the stiffness /3. The equivalence of the non-linear o model with the 
lagrangian in Eq. © implies that the stiffness evaluated in the two cases must coincide. 
Therefore one can insert the expression of /3 of Eq. Q, properly modified to account for 
the finite T and /! effects, i.e. 



P = 2g 2 N f N c T £ / 

n — — oo <> A. 
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into Eq. © and determine the corresponding solution for a or A. 

As an example, in fig. 2 we plot the solutions of Eq. ©: a and the mass of the pionic 
mode m n = \fk versus the temperature, at jl = and for two different values of the 
number of colors: N c = 3 (solid lines) and N c = 10 (dashed lines). The values of the two 
cutoffs entering Eqs. © and (fTOb are A n = 200 MeV and A = 675 MeV. In fact A^ and 
A are related to different degrees of freedom and do not need to have the same value and 
we take a value of A n close to f n , which fixes the scale of the non-linear a model. (The 
choice A^ ~ A made in has been criticized in [4]. For other references on this point 
see Ja].) 

We observe in fig. 2 that at low temperatures the solution of Eq. © corresponding 
to A = (not plotted) and non-vanishing a is realized. Then o decreases with T and 
eventually vanishes. By further increasing the temperature, the solution o = (not 
plotted) and m n = ^ is realized. 

The two regions meet at T = T c where o = A = 0. This corresponds in Eq. © to the 
critical stiffness 
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and T c is determined by the comparison of Eqs. (flOb and (ITTb In the particular case 
considered in fig. 2 we get T c w 146 MeV for N c = 3 and T c w 172 MeV for N c = 10. 

One should notice that in the non-linear o model the region below T c shows a non- 
vanishing expectation value of the field o which implies a dynamical breaking of the 
chiral symmetry. Above T c this expectation value vanishes and at the same time a finite 
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FIGURE 2. a and m n vs. T at = (in MeV) as obtained from Eq. ©. See text for the specific values 
of the other parameters. 



mass for the pions is generated. The absence from the spectrum of massless Goldstone 
bosons is a signal of chiral symmetry restoration. Therefore, within this framework, T c 
is the critical temperature related to the chiral symmetry. 

In fig. 2, when the temperature is increased, m n grows and, at some temperature, 
Tpair, eventually diverges and the % mesons decouple and disappear from the spectrum. 
Reasonably, T pc dr could be interpreted as the temperature of dissociation of the qq pair. 
A remarkable feature which appears in fig. 2 is that T va x r does not depend on N c and, 
in practice, it numerically coincides with the critical temperature T* determined by the 
mean field analysis: Tpair ~ T* ~ 185 MeV. 

Therefore, the comparison of the NJL model with the non-linear o model indicates 
that fluctuations indeed modify the mean field conclusions and the new critical temper- 
ature T c , associated with the chiral symmetry restoration, is generated. T c turns out to 
be smaller than T* so that we find a temperature interval T c < T < T* where the gap 
equation still provides m q ^ but the chiral symmetry is restored as an effect of the 
fluctuations. Only above T* ~ Tpair the mesons disappear and m q vanishes. Consistently 
with this picture we note that T c depends on N c and T c — > T* when N c — > °°, i.e. when 
the fluctuation effects become suppressed and the mean field regime holds. 

A comment about m q is in order. In fact the common belief is that a non-vanishing 
m q is directly related to the breakdown of chiral symmetry. This is true only if the 
gap m q appears as a pole in the two point function of the fermion i//. In our specific 
case m q corresponds to the pole mass of ty/ only in the mean field approximation. This 
property has been first observed in 2D models in fl2ll . and, for a detailed treatment of 
the problem, we refer to 11131 \\M . With a suitable change of variables, corresponding 
to the introduction of polar 'coordinates': o + iy 5 ^ = p exp(/y 5 0) and of the fermion 
X = exp(?'y 5 0/2), where p and x are chirally neutral, it is possible to show that m q , 
as obtained from the gap equations and © is a pole in the x two point function 
(of course in our specific case the definition of the new variables should be modified 



to account for the flavor degrees of freedom). Due to the chiral neutrality of %, this 
excitation does not break chiral symmetry. Conversely, it does not appear in the \j/ 
propagator as a pole but rather, due to the 6 field fluctuations, as a branch cut. In this 
sense, a finite m q is still compatible with chiral symmetry. 

The presence of a region between two characteristic temperatures where the fluctua- 
tions cancel the long range order associated to the quark mass, restoring the symmetry, 
shows a clear resemblance with the pseudogap phase of high temperature superconduc- 
tors 10. 0,113] and this motivates the following analysis of the spectral function. 



EVALUATION OF THE PSEUDOGAP 

After having determined the two characteristic temperatures T c and T*, we conclude 
by showing that a simple approximation in the calculation of the spectral function 
provides clear indications of a pseudogap Ja]. The spectral function of the fermionic 
quasiparticles, N(co), is obtained from the imaginary part of G R (k, ft)): 

N{(0) = -\S (2^ Tr ^ rf ^ Im GR{% > W) (12) 

where Tr c ,f,d is the trace over Dirac, color and flavor degrees of freedom. G R (k, ft)) is 
the analytical continuation of the imaginary time retarded Green function: G R (k, ft)) = 



G (k , ft) + z'e) — ft)) where G (k, ft) + z'e) corresponds to the free term con- 
tribution and ft)) includes higher order corrections. The starting point to calculate 
the retarded Green function is the Lagrangian 

^ eff = ^{id v f + m - gQ {o + iy 5 x-K)}w-^[o 2 + Ti 2 ] (13) 

where the o field is replaced by expanding the constraint : o = (f% — 7t 2 ) 1//2 ~ f n — 

7i 2 /(2f n ). The correction £/?(fc, ft)) to the free fermion propagator is computed to the 
one loop level from Eq. (fT3l . [11], and its imaginary part generates the correction, N pert , 
to the free fermion contribution in the spectral function 

N((0) = ^f{(D + ll)J{(Q + ll) 1 -m 2 q + Npert ■ (14) 

In fact we have included the mass m q , which is dynamically generated, in the free 
fermion contribution and we have also included in the definition of go the form factor 

So(<l 2 ) = So m p/ (mp+ (2nnT) 2 + q 2 ) 2 with mp = 100 MeV, which acts as a regulator 
that corrects the large momentum behavior of the loop. 

We normalize N(co) to the free massless case : R(co) = N{(0) / {^^-{(0 + /J,) 2 ) and 
plot R(co) for (T-T c )/T c = 0.02,0.15,0.20 (ju = 10 MeV) in fig. 3. The solid straight 
line corresponds to T = T*. In fig. 3 we observe the pseudogap for values of ft) just above 
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FIGURE 3. The ratio N((o)/Nf ree {(o) as a function of the energy (0, at /! = 10 MeV and for various 
temperatures between T c and T* . See text. 



the value CO = m q — /i, where the free part of N(co) vanishes, whereas it disappears for 
large values of CO. The effect of the pseudogap is larger for temperatures just greater 
than T c and vanishes at T = T*. Therefore, even in this simple approximation, it is 
possible to realize the presence of the depletion of states in the spectral function and this 
calculation supports the picture of a pseudogap appearing between the two characteristic 
temperatures. Since the NJL model effectively contains many ingredients of QCD, this 
analysis is a starting point to check the existence of a pseudogap in QCD and its possible 
role in the explanation of the mesonic bound states observed in lattice calculations. 
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